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We study systems of multiple interacting quantum impurities deposited on a metallic surface in 
a 3d host. For the real-space two-impurity problem, using numerical renormalization group calcu¬ 
lations, a rich range of behavior is shown to arise due to the interplay between Kondo physics and 
effective RKKY interactions — provided the impurity separation is small. Such calculations allow 
identification of the minimum impurity separation required for a description in terms of indepen¬ 
dent impurities, and thereby the onset of the ‘dilute impurity limit’ in many-impurity systems. A 
‘dilute cluster’ limit is also identified in systems with higher impurity density, where inter-impurity 
interactions are only important within independent clusters. We calculate the quasiparticle inter¬ 
ference (QPI) due to two and many impurities, and explore the consequences of the independent 
impurity/cluster paradigms. Our results provide a framework to investigate the effects of disorder 
due to interacting impurities at experimentally relevant surface coverages. 
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I. INTRODUCTION 

Electronic scattering due to impurities in metals gives 
rise to a variety of experimental signatures.^ Bulk mea¬ 
surements, such as resistivity, are of course affected; but 
local quantities provide much more information about 
the effect of impurities. In particular, the development of 
scanning tunneling spectroscopy (STS) has allowed im¬ 
purity positions/distributions on the surface to be im¬ 
aged directly, and real-space maps of the surface local 
density of states (LDOS) to be built up.^ These maps 
show pronounced inhomogeneities and even quantum in¬ 
terference due to the presence of impurities.^ 

More recently, ET-STS - the Eourier transform of 
LDOS maps obtained through STS - has been connected 
to quasiparticle interference (QPI).^ It describes how the 
Bloch state quasiparticles of a clean host system are scat¬ 
tered from impurities. Experimental QPI patterns are 
typically interpreted in terms of scattering from a single 
potential defect in the weak-scattering (Born) limit. 
The clean host structure can then be deduced from the 
preferred QPI scattering vectors. 

The phenomenon of impurity-induced scattering is 
much richer in the case of magnetic impurities. De¬ 
scribed in terms of a single interacting quantum level, 
such an impurity can host a free local moment at high 
energies/temperatures, which is then screened in metal¬ 
lic systems by conduction electrons through the Kondo 
effectbelow a characteristic scale The en¬ 

hanced spin-flip scattering responsible for the Kondo ef¬ 
fect at low energies leads to a characteristic impurity 
spectral resonance^observable in STS.^ Induced 
real-space LDOS modulations from a single magnetic im¬ 
purity produce a strong QPI signal, which acquires a 
universal temperature and scanning-energy dependence 
at low energiesdue to universality of the Kondo res¬ 
onance in terms of and In the STS 


context, the scanning-energy uj = eVsd is simply related 
to the bias voltage Vsd between STM tip and surface.^ 
The dynamics of the QPI therefore provides information 
on the scattering t-matrix of the system. 

Unlike the case of a simple potential defect, interacting 
quantum impurities give rise to non-trivial many-body 
effects. The numerical renormalization group^’^^ (NRG) 
is the theoretical method of choice for treating a single 
quantum impurity: it provides numerically-exact access 
to dynamical impurity quantities on essentially any tem¬ 
perature/energy scale.Solution of the local impurity 
problem then allows other real-space quantities'^ and the 
QPI to be calculated, as detailed generally in Ref. 12. 

However, real systems contain many impurities: the 
real-space surface region of a host material probed by 
STS to obtain QPI could include N ^ 10—100 impurities, 
depending on sampling size and impurity density. Eor 
potential defects, the physics remains rather trivial (the 
scattering t-matrix can be obtained in closed form^^). 
By contrast, systems containing many magnetic impu¬ 
rities embody subtle interplays between Kondo physics, 
effective through-host RKKY interactions, and disorder 
in the impurity distribution. This is the scenario studied 
in the present work. 

We specify a model in Sec. II that takes into ac¬ 
count explicitly the surface of a 3d metallic host onto 
which multiple magnetic impurities are deposited, allow¬ 
ing thereby calculation of surface LDOS maps and QPI. 
Realistic parameters are chosen, yielding single-impurity 
Kondo temperatures in the range ^ q 1 _ lOOK rel¬ 

evant to experiment. Electronic and quasiparticle scat¬ 
tering in this system is described in terms of the t-matrix, 
itself controlled by impurity dynamical quantities. The 
theoretical prescription for calculating QPI for a many- 
impurity system in terms of the impurity Green functions 
is given in Sec. Ill, including (exactly) the effect of sur¬ 
face quasiparticle dephasing by the host bulk. 
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Solution of such a model with N interacting quantum 
impurities is a formidable challenge because in general 
it involves N coupled screening channels. Generalization 
of NRG to deal with two spatially-separated impurities 
already represents a very significant increase in computa¬ 
tional complexity, although accurate results for dynam¬ 
ical quantities are now possible. Our approach to the 
many-impurity problem in this paper begins with a de¬ 
tailed analysis of the two-impurity model, solved using 
NRG in Sec. IV. Inter-impurity interactions are impor¬ 
tant when the impurities are separated by a few lattice 
sites, giving rise to a rich range of correlated electron 
physics with distinctive QPI signatures. 

A key finding, however, is that the impurities behave 
essentially independently for surprisingly small inter¬ 
impurity separation. This has implications for the on¬ 
set of the ‘dilute limit’ in many-impurity systems. In 
this case, local quantities are in essence ‘blind’ to the 
impurity distribution (although the impurity structure 
factor does affect the QPI). By analyzing random impu¬ 
rity distributions on the surface in Sec. V, we find that at 
experimentally-relevant impurity concentrations up to a 
few %, there is a scale separation in the occurrence of im¬ 
purity clusters of a given size Nc. This allows a descrip¬ 
tion of many-impurity systems in terms of independent 
impurities and independent elusters. As a representa¬ 
tive example, we consider V = 100 impurities randomly 
distributed in a 500 x 500 surface sample region of a semi¬ 
infinite 3d cubic lattice host. On average, three Nc = 2 
clusters appear in such a sample: they are well-separated 
from other impurities or clusters and can be treated in¬ 
dependently; and no clusters of size Nc > 2 appear on 
average. Intra-cluster interactions are important for cap¬ 
turing local spectroscopic details; but the overall effect 
on QPI is weak due to low cluster occurrence at typical 
impurity densities. 

II. MODEL FOR MULTIPLE MAGNETIC 
IMPURITIES 

We consider N > 1 magnetic impurities deposited on 
the surface of a metallic host system. The full Hamilto¬ 
nian is given by 

N 

H — -Hhost T ^ ^ -f^impjCk: • (l) 

OL = l 

The host is taken here to be an infinite 3d cubic tight 
binding lattice, cleaved to reveal the (100) surface, 

CX) 

i^host = -* V [ V + H.C.) 

Z=0 {ij),a ^ 2 ) 

^ 5^(^ri2;cr^ri(2+l)cr H.C.) , 

i,a 

where creates an electron of spin a =tA in the sin¬ 
gle Wannier orbital localized at site of layer z, and (ij) 


denotes the sum over nearest-neighbor sites of a 2d square 
lattice layer. The resulting LDOS on the translationally- 
invariant surface is finite and fiat at low energies, so al¬ 
though we focus on this specific host realization, the real- 
space physics is rather generic and typical of metallic 3d 
systems. More complicated materials can always be cast 
in the layered form of Eq. 2 by generalizing to a matrix 
structure. 

We treat the surface explicitly to allow study of the 
QPI, which is necessarily obtained through surface mea¬ 
surement in FT-STS experiments. We note that coupling 
the surface layer to the semi-infinite bulk is required to 
capture the important effects of surface quasiparticle de¬ 
phasing. Properties of the clean host are characterized 
by its real-space Green functions, here calculated exactly 
using the convolution method introduced in Ref. 12. 

Each magnetic impurity a is described in terms of a 
single interacting quantum level, tunnel-coupled to a sur¬ 
face {z = 0) host site at position 

+ H 22 + H.C. j , 

a 

where creates a spin-a electron on impurity a. Im¬ 
portantly, such an Anderson impurity has internal spin 
and charge dynamics due to the local electronic inter¬ 
actions. This is in contrast to a static ‘magnetic’ im¬ 
purity caricatured by an inhomogeneous magnetic field 
that breaks underlying time-reversal and spin SU(2) sym¬ 
metry (and is non-interacting). The generalization of 
Eq. 3 to multi-orbital impurities is not considered in the 
present work (the formal aspects of the scattering prob¬ 
lem being largely unchanged). 

The generalized quantum impurity problem involving 
N magnetic impurities, spatially separated and coupled 
to conduction electrons of the host lattice, is distinctly 
nontrivial and exhibits a rich range of correlated electron 
physics. The simplest N = 1 model, comprising a sin¬ 
gle magnetic impurity, already exhibits strong dynamical 
effects due to interactions, such as the Kondo effect. 

In the case of multiple impurities, richer physical be¬ 
havior arises due to an effective RKKY interaction which 
couples impurities indirectly. The RKKY interaction 
is mediated via the host lattice and therefore depends 
on the specific impurity distribution in real-space. As 
shown explicitly in Sec. IV B, impurities can be Kondo 
screened independently by the host, or entangled clus¬ 
ters can be collectively screened in a multi-stage process. 
Impurities can also ‘screen themselves’ by forming inter¬ 
impurity singlet states when the effective RKKY inter¬ 
action is strongly antiferromagnetic. 

A. Impurity parameters 

The Kondo physics of a single magnetic impurity is 
sensitive to the underlying model parameters and host 
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material through its LDOS.^’^^ However, the RG flow 
and associated universality is controlled by an emergent 
energy scale — the single-impurity Kondo tem¬ 
perature. In the Kondo limit,^ determined by 

the interaction strength U/Tq^ where Fq = and 

po = l/6t is the Fermi level surface LDOS of the 3d cu¬ 
bic lattice. 

In real systems, Kondo temperatures measured in ex¬ 
periment are known to vary widely, even for a given im¬ 
purity type and host material.This is principally due 
to differences in the hybridization Fq, which is sensitive 
to details of the impurity’s local environment. 

Guided by this, in the present work we use realis¬ 
tic, fixed values of the conduction electron bandwidth, 
12t = 11 eV, and impurity interaction strength U = 
3t = 2.75 eV (taken from Refs. 22 for Go impurities on 
Au). For convenience, we consider Cd = —U/2 such that 
the impurity is singly-occupied. Appropriate choice of 
V then gives rise to a realistic spread of single-impurity 
Kondo temperatures, as summarized in Table I. 


TABLE I. Kondo temperatures for 
a single impurity on the 3d cu¬ 
bic lattice surface, with bandwidth 
12t = 11 eV and U = 2.75 eV. Cal¬ 
culated as the half-width at half¬ 
maximum of the Kondo spectral 
resonance at T = 0 using NRG. 


We note that U/Tq in the range 12-16 yields Kondo 
temperatures consistent with classic studies of Go atoms 
on a Gu or Au surface;^’^^’^^”^^ while a passivating layer^^ 
of Gu 2 N between surface impurities and bulk reduces the 
hybridization to yield ^ 2K. Indeed, U/Tq = 24- 

28 is more appropriate for Fe in Au, where ^ 

0.3K.2^ 

In the following we use these impurity parameters in 
the context of many-impurity systems, where different 
emergent energy scales and physics naturally arise. 


f//ro 

ry-ilimp 

12 

85K 

16 

15K 

20 

3K 

24 

0.5K 

28 
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III. ELECTRON AND QUASIPARTICLE 
SCATTERING 

The impurity single-particle dynamics are described 
generically by elements of the Green function matrix 
[Gd{ui)]a ,/3 = G'fiu:) = where {{A;B))^ 

is the Fourier transform of the retarded correlator 
—i0{t){{A{t),B{0)}). The propagator therefore 

contains information on the energy-dependent scattering 
of electrons between impurities a and /3. 

For a many-impurity system, such Green functions are 
obtained from a matrix Dyson equation, 

= (4) 

where the non-interacting (but host-coupled) impurity 
Green functions are given by 

[gd(w)]“^ = (w + iO+ -ed)I-r(w) , (5) 


in terms of the hybridization matrix T{uj) with elements 
[r(w)]a ,/3 = V^GlQiYocYjd.uj). Here, G°^,(r„,r; 3 ,w) = 
((Cr„ 2 <T; 4 ; 32 v))° is the propagator between site r„ of 
layer z and of layer z' in the clean host (without impu¬ 
rities). and are the host sites to which impurities 
a and {3 are coupled on the surface layer, z = 0. The 
self-energy matrix 5](co’) contains all information due to 
electronic interactions, which give rise to the Kondo ef¬ 
fect, RKKY interaction, etc. 

Likewise, the full Green functions de¬ 

scribe electronic propagation through the host, in the 
presence of impurities. The contribution to the total 
scattering from the impurities is given by the real-space 
t-matrix equation, 

G^^,{Yi,Yj,Uj) =Gl^,{Yi,Yj,Uj) 

+ Ya,Uj)Tal3{uj)Gl^, {y^, Yj,Uj) 

Ck:,/3 

where the sum runs over impurity sites a and jS; and 
the t-matrix is given explicitly in terms of the impurity 
Green functions by 

T„^iuj)=V^Gfiuj). (7) 

While host surface sites are of course coupled to those 
in the bulk, the STM tip itself probes only the host sur¬ 
face. In consequence, as detailed in Ref. 12, the scat¬ 
tering problem is most effectively formulated in the so- 
called surface diagonal basis, in which the layer index is 
preserved; leading to a description of the QPI in terms of 
surface quasiparticles. This follows from the partial diag- 
onalization of iFhost (Eq. 2) by 2d Fourier transformation 
of layers parallel to the surface. 


^kii za 


i-k|| 


Q 


1/2 

BZ Vi 


(8) 


where is the volume of the first (surface) Brillouin 
zone (IBZ). The host Hamiltonian then reduces to a bun¬ 
dle of decoupled Id chains. 


-f^host ^ ^ ^ ^ ^k|| ^k|| 2:0-^k|| za ^(^ky 2;cr^k|| (2;+l)crd'E-C.) , 

k||,cr z =0 

(9) 

where = —2t[cos{aQkx) + cos(aoA:y)] is the 2d square 
lattice dispersion, given in terms of surface momentum 
k|| = {kx^ky), and with lattice constant oq. As a conse- 

quence, G°(k||,k||,w) = <^(k|| - kj|)- 

Surface quasiparticles are however dephased by the bulk 
even in the clean host,^^ with 


GUdh,L0) = f 


UJ — €\^2d 

^11 

2t 


where 


= UJ 


|^sgn(w)-\/w^ — 1 


( 10 ) 


Wl — GT 


|w| > 1 
|w| < 1 
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Scattering of surface quasiparticles due to impurities is 
then described^^ by a surface t-matrix equation, 


Go(k||,k||,a;) =G0(k||,c^)5(k||-k') 

+ G°(k||,a;)T(k||,kj|,c^)G'](k||,a;), 

( 11 ) 


with the t-matrix itself given by, 

r(k||,kj|,:^) = ^ , ( 12 ) 

^ ^BZ Q 

a,l3 

A. Quasiparticle interference 


Quasiparticle scattering from impurities gives rise to 
QPI, as measured in experiment by FT-STS.^ It is ob¬ 
tained from the LDOS map p(r^, cj) at scanning-energy uj 
and temperature T, measured in real-space over a sample 
region of size L x 

P(q>w)= XI . (13) 

ie (^L X Z/) 

For magnetic impurities deposited on the (100) surface 
of a semi-infinite 3d cubic lattice, the QPI can be calcu¬ 
lated following Ref. 12. Taking the sample size L ^ oo 
and subtracting the trivial contribution at q = 0 from 
the clean system, the exact QPI, Ap(q, cj) = p(q, cj) — 
p^(q, cj) from Eq. 13 is given by, 

Ap(q, ^) = - ^ [e(q> - Q(-q. w)*] > (14) 

where 

Q(q,w)= [ [Go(k||,k|| - q,w) - G^(k||,w)<5(q)] 

J ^‘‘BZ 

IBZ (15) 

= V^Y^Gf{w)xKp{ci,oj) 

a ,13 

is given in terms of the full impurity Green functions 
The quantity Aci/ 3 (q, cj) is defined for the clean 
host material, and follows from Eqs. 11 and 12 as, 

Ao^{q,w)= f Go(k||,w)Go(k||-q,w) 

dz ( 1 «) 

which can be computed accurately and efficiently as a 
convolution. 


IV. TWO IMPURITIES 

The two-impurity problem is the simplest to cap¬ 
ture the competition between Kondo physics and inter¬ 
impurity interactions.A leading-order perturbative 


treatment generates the host-mediated RKKY interac¬ 
tion, coupling impurities by indirect exchange. 

The magnitude and sign of the RKKY interaction de¬ 
pends sensitively on the dimensionality and geometry 
of the host,^^ and on the impurity separation vector 
R = r2-ri.30 

Two spatially-separated Anderson impurities are often 
modelled by a two-impurity Kondo model ( 2 IKM),^^’^^ 
describing two exchange-coupled spin-| impurities each 
coupled to its own independent conduction electron chan¬ 
nel. The physics of this model is immensely rich: Kondo- 
screened and inter-impurity singlet phases arise, sepa¬ 
rated by a non-Eermi liquid quantum critical point 
of two-channel Kondo type.^^ 

The 2IKM is however oversimplified, because the two 
conduction electron channels are not strictly indepen¬ 
dent in a true real-space system, both being constructed 
from states of the same electronic host.^^ This is reflected 
by finite off-diagonal elements of the hybridization ma¬ 
trix, T{uj) in Eq. 5. Since ri 2 (co’) = K^Goo(i*i, r 2 , cj) 
involves the host propagator between impurities located 
at sites ri 7 ^ r 2 , the two-impurity problem evidently 
features inter-channel charge transfer processes. These 
processes are RG relevant: simple perturbations that 
break the same symmetries can be incorporated directly 
into the 2IKM, and are known to destabilize the critical 
point.Low-temperature/energy Eermi liquid crossovers 
then arise ubiquitously.^^ 

The true real-space model does not therefore support 
a quantum phase transition. Instead, there is a crossover 
as a function of impurity separation or impurity-host 
coupling, with the Kondo regime evolving continuously 
into an RKKY-dominated regime.This crossover 
is also known from the simpler two-impurity Anderson 
model, in which a hopping t' directly tunnel-couples the 
impurities.In that case, the RKKY interaction is 
generated to second-order in the tunnel-couplings, while 
the inter-channel charge transfer arises to third-order.^^ 

It is important to emphasize that inter-channel charge 
transfer processes cannot simply be neglected: finite 
Goo(i*i, r 2 , cj) is the common origin of both the RKKY 
interaction itself and the relevant perturbations destroy¬ 
ing 2IKM criticality. Eurthermore, the dynamical na¬ 
ture of Goo(i*i, r 2 , cj), resulting from through-lattice elec¬ 
tronic propagation, produces nontrivial RG flow. Phys¬ 
ical regimes of the real-space model might therefore be 
inaccessible within the 2 IAM (where Goo(ri, r 2 , cj) is re¬ 
placed by a real constant hopping P), or in the 2IKM 
(where it is neglected altogether). 

As such, if one is interested in the real-space physics 
of two impurity systems — and the resulting QPI — one 
must study a real-space model directly. Here we solve 
the true real-space two-impurity problem exactly using 
NRG,^’^^ as described below. 
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A. Numerical Renormalization Group 

For two impurities separated in real space on the 
host surface, the first step is to diagonalize the matrix 
Dyson equation, Eq. 4. The translational invariance of 
the host surface implies that this can be achieved for 
all cj by a single canonical transformation of operators 
to an even/odd orbital basis. This yields = 

{uj + i 0 + — ed — re/o(co’) — in terms of even/odd 

quantities fte/o = i ^12 (noting that Qn = Q 22 
and Cli 2 = In the non-interacting system {U = 0) 

where T^e/o = 0 , the even and odd channels are strictly 
decoupled: the even(odd) impurity combination couples 
only to the even(odd) host combination. However, the 
even/odd transformation scrambles the interaction term 
in coupling even and odd channels for U 7 ^ 0, 

thereby requiring an irreducible two-channel NRG cal¬ 
culation. The original information about the real-space 
separation of the impurities is encoded in the differ¬ 
ence between the dynamical quantities re(co’) and ro(co’). 
As noted in Refs. 33, 30, the full energy dependence of 
these functions is required to capture the true real-space 
physics. 

The explicit transformation of real-space conduction 
electron operators to an even/odd basis (as described in 
Ref. 28) can be highly complicated, depending on the 
host lattice and inter-impurity separation vector. We 
note however that this is not required, since the transfor¬ 
mation can be performed on the level of the continuous, 
energy-dependent hybridization functions, as above. 

The two-impurity NRG calculation then involves a 
logarithmic discretization^’of even/odd hybridization 
functions A discretized version of the full model 

is then formulated in terms of even and odd impurity 
combinations, coupled to the end of even and odd Wilson 
chains. The model is then diagonalized iteratively, start¬ 
ing from a subsystem comprising the impurities them¬ 
selves, and then building up the chains by successively 
coupling on even and odd Wilson chain orbitals. The 
couplings down each Wilson chain decrease exponentially 
due to the logarithmic discretization, and so high-energy 
states can be safely discarded at each step to avoid expo¬ 
nential Hilbert space growth. The physics of the problem 
is revealed on progressively lower energies/temperatures 
as the Wilson chains grow in length — this is the essential 
RG character of the quantum impurity problem. 

The computational complexity of real-space two impu¬ 
rity problems is significantly greater^^ than that of single 
impurity problems, which involve only a single effective 
conduction electron channel; although accurate results 
are now possible. Further, impurity self-energies can be 
calculated accurately within NRG through the identity, 

S(w) = [Gd{oj)]-^Fd{oj) , (17) 

where [Fci{ui)]a ,/3 = Up-))^ (a matrix gen- 
eralization to many impurities of the result originally 
given for the single-impurity Anderson model^^). Both 


Gd{w}) and Fc^(co’) are calculated directly in NRG using 
the full density matrix approach^ within the complete 
Anders-Schiller basis.Impurity Green functions then 
follow from the matrix Dyson equation, Eq. 4. 

In the following, the real-space two-impurity prob¬ 
lem is solved with NRG using a discretization param¬ 
eter A = 2 , and retaining Nk = 15,000 states at each 
step of the iterative process. Gonserved total charge and 
spin projection were exploited for block diagonalization. 
Gonvergence was established by increasing A and Nk- 
^-averaging was not required. 


B. Effect of impurity separation 


The interplay between Kondo physics and the effective 
host-mediated interaction between two magnetic impuri¬ 
ties on a surface, gives rise to a rich range of physics. 
Aside from the bare energy scales t, U and V enter¬ 
ing the problem, two emergent scales arise: the single¬ 
impurity Kondo temperature and the effective 

inter-impurity RKKY exchange coupling, Jrkky- The 
relative size of these scales, and the sign of Jrkky, con¬ 
trols the underlying RG flow and low-energy physics. 

In the following, we use NRG to solve the two impurity 
problem, allowing these scales to be extracted for given 
impurity parameters and inter-impurity separation vec¬ 
tor R. However, it is also instructive to obtain simple 
perturbative estimates. For a single impurity on the cu¬ 
bic lattice surface, third-order perturbative scaling^at 
particle-hole symmetry yields. 



where Tq = 7 rE^/ 6 t, as in Sec. HA. The estimates com¬ 
pare well with exact (non-perturbative) results in Table I 
using a ^ 3. For two impurities, the RKKY coupling is 
given to second order vhTq/U as,^^ 




where P denotes the principal value, and Goo(R, c) = 
Goo(ri,r 2 ,e) are free host Green functions connecting 
surface impurity sites ri and r 2 . At large inter-impurity 
lattice separations, we find from Eq. 19 that the RKKY 
coupling decays rapidly, faster than l/|Rp. 

As such, when the impurities are well-separated, one 
expects ^ Jrkky for a given U/Tq. Each impu¬ 

rity is then individually Kondo-screened essentially in¬ 
dependently of the other, thereby quenching their mag¬ 
netic moments and rendering the RKKY exchange inter¬ 
action inoperative. By contrast, at small inter-impurity 
separations such that Jrkky ^ RKKY physics 

should dominate, causing the impurities to bind together 
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FIG. 1. Two equivalent mag¬ 
netic impurities on the 3d cu¬ 
bic (100) surface with U/Fq = 
28 (see Table I). NRG results 
for the local impurity spec¬ 
trum —lm[tGy{(jj)] vs (jj/t at 
T = 0 (upper panels) and the 
temperature-dependence of the 
impurity contribution to en¬ 
tropy 5'imp(T) vs T/t for the 
same systems (lower panels). 
Inter-impurity vector R, with 
= 0 (left panels) and R^ = 
Ra; (right), where Ra;/ao = 
1, 2, 3,4 for black, red, blue and 
green lines. Gircle points are 
the single-impurity result. 


in either singlet or triplet configurations (depending on 
whether Jrkky is positive or negative). In this scenario, 
the single-impurity RG flow is cut off on the scale of 
-^RKKY, and single-impurity Kondo physics is destroyed. 
We remark however that, even in the antiferromagnetic 
case, the inter-impurity singlet remains entangled with 
the host, which mediates the indirect exchange coupling 
(in contrast to the situation in the 2IKM with direct ex¬ 
change, where the singlet decouples). 

TABLE II. 

Perturbative estimate 
of the RKKY coupling 
from Eq. 19, given as 
(e/RKKv/t) X (fZ/To)^, 
for inter-impurity sepa¬ 
rations R = (Rx^Ry) as 
per Eig. 1. 

Table II shows the perturbative estimate for the 
RKKY coupling, as obtained from numerical evaluation 
of Eq. 19, for various inter-impurity separations. For sep¬ 
arations with Ry = 0, there is a clear alternation between 
antiferromagnetic and ferromagnetic coupling; while for 
Ry = Rx, only ferromagnetic coupling arises. The rapid 
decay of Jrkky with increasing separation |R| is also 
clearly seen (even for small separations) from this sim¬ 
ple calculation. The true variation of Jrkky — and its 
subtle interplay with Kondo physics — must however be 
extracted from the full NRG solution of the two-impurity 
model, which takes into account the dynamical nature of 
the coupling and renormalization effects. We now turn 
to NRG to examine these issues in detail. 

Physical quantities, such as impurity dynamics and 
thermodynamics, provide detailed information about 
the underlying physics and screening mechanisms. In 
Fig. 1 we use NRG to study the impurity spectrum 
—Im[tGy (cj)], plotted vs uj/t at temperature T = 0 
in the upper panels; and the impurity contribution to 
entropy 5'imp(J")/te, shown vs temperature T/t in the 


Rx 

II 

o 

II 

1 

+4.5 X 10"^ 

-2.5 X lO-"" 

2 

-2.5 X 10-3 

-2.0 X 10-3 

3 

+1.0 X 10"^ 

-2.9 X 10-3 

4 

-2.5 X 10"® 

-3.2 X 10-3 


lower panels. The two impurities are deposited on the 
cubic lattice surface, separated by a vector R = (R^,, R^); 
with R^ = 0 considered in the left panels and R^ = R^, 
in the right panels, and with R^^/ao = 1,2, 3,4 in each 
case. Note that GY{uj) = since the impurities 

are equivalent. The local spectra are related to the local 
parts of the scattering t-matrix through Eq. 7. 

We consider first the spectra in the upper panels. 
In all cases, a Hubbard satellite peak arises at high 
energies |ci;| ^ G/2, corresponding to impurity charge 
fluctuations.^ For |cc;| ^ 17, the impurities become essen¬ 
tially singly-occupied and magnetic moments form. Al¬ 
though the impurities are coupled by the effective RKKY 
interaction, at energies or temperatures ^ Jrkky, 
the impurities behave essentially independently. Cor¬ 
rections to the local moment fixed point then control 
the physics; incipient RG flow for each impurity to¬ 
wards a Kondo-screened strong-coupling state produces 
the classic inverse logarithmic spectral behavior^ when 

max(| JrkkyI, ^ ^ U, 


-lm[tGy{uj)] 


1 

lTb\n^\cuj/T^^^\ ’ 


( 20 ) 


with h^c constants of 0(1). 

Consider now the R^ = 0 cases (left panels. Fig. 1). 
For neighboring impurities (R^^/ao = 1, black lines) the 
RKKY interaction Jrkky > 0 is antiferromagnetic and 
comparatively strong, causing the impurities to lock up 
together into a spin singlet state on the scale of Jrkky- 
The RG flow toward independent Kondo strong coupling 
states is arrested on this scale since all impurity degrees 
of freedom are quenched. This results in an impurity 
spectrum lm[tGY] — 0 for jcjl ^ Jrkky- Confirming 
this physical picture, the impurity contribution to en¬ 
tropy shown in the lower panel is quenched to 5'imp = 0 
on the same scale (for T ^ | Jrkky |, S'imp ^ ^B ln(4), 
corresponding to two free impurity local moments). The 
effective RKKY coupling in this system can thus be ex¬ 
tracted from these crossovers as Jrkky ~ 
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By contrast, when R^c/ao = 2 (left panels, red lines), 
the effective RKKY interaction is ferromagnetic and 
weaker. On the scale of |Jrkky| ~ the impuri¬ 

ties form a triplet configuration with an associated ln(3) 
impurity entropy. On lower temperature/energy scales 
T|=i 10-^t, we find that the inter-impurity triplet 
is exactly Kondo screened^^ by even and odd conduc¬ 
tion electron channels. On energy scales Ic^l ^ 

|-^rkky|, spin-1 local moment physics naturally arises, 
with spectral behavior again that of Eq. 20, but with the 
characteristic scale replaced by In fact, uni¬ 

versal lineshapes are obtained as a function of \uj\/T^^^ 
and T for all |cc;| and T <C | JrkkyI- In particular, 
the lowest-energy behavior is characteristic of flow to a 
Fermi liquid ground state, with a fully quenched entropy, 
5'imp(0) = 0. The low-energy behavior of the impurity 
spectrum has the characteristic Fermi liquid asymptotics, 

-MtGYiuj)] ’ ( 21 ) 

where Tk = is the triplet Kondo scale, q = 0(1) 

is a constant, and the Fermi level value p depends on 
the full self-energy of the two-impurity problem (there is 
no simple pinning condition, as for the single-impurity 
case).^^ 

On further separating the impurities, the RKKY inter¬ 
action becomes even weaker. In fact, |Jrkky| ^ ~ 

10“^t already for R^^/ao = 3, 4 - left panels, blue and 
green lines (which are virtually indistinguishable) - such 
that each impurity is independently Kondo screened be¬ 
low the single-channel scale The lineshapes are 

again universal,^ although different from the triplet case 
in detail. Comparison with the circle points shows that 
the physics of each impurity is indeed that of a single im¬ 
purity in the same host. Here the impurity spectrum is 
also described by Eq. 21, but now with p = To/t pinned 
by the Friedel sum rule.^’^^ 

Along R^ = Ra, (Fig. 1, right panels), one might expect 
the effective RKKY coupling to be ferromagnetic, since 
the impurities are always connected by an even number 
of hops on the lattice. As demonstrated by Fig. 1, this is 
indeed the case for R^c/ao = 1 and 2. An inter-impurity 
triplet thus forms in both cases, with incipient RG flow 
near 5'imp(T) = /cBln(3). For R^c/ao = 1 (black lines), 
the Kondo screening of the inter-impurity triplet is 
distinctly two-stage^ indicating a strong inequivalence of 
the even and odd conduction electron channels. This 
results from significant off-diagonal elements of the hy¬ 
bridization matrix r(u;) in Eq. 5 when the impurities are 
proximal. However, this anisotropy is already very small 
when Ra;/ao = 2 (red lines), so that the two-channel 
Kondo screening of the inter-impurity triplet proceeds 
essentially in a single step. At greater separations 
Rcc/ao = 3 and 4 (blue and green lines), the magnitude 
of the RKKY interaction is again small enough that the 
two impurities remain essentially independent, and are 
separately Kondo screened. 



FIG. 2. Evolution of impurity spectra —roImGj^(ct;) with de¬ 
creasing H/Fo, obtained by NRG. Plotted at T = 0 vs uj/t for 
Rcc = R^ = 1, 2, 3,4ao as the black, red, blue and green lines. 
Gircle points are the corresponding single-impurity results. 


For two impurities on the surface of a 3d cubic lat¬ 
tice, the overall results of Fig. 1 shows that the dy¬ 
namics are strongly affected when their separation is 
on the order of a few lattice sites. Importantly, how¬ 
ever, the impurities are found to behave essentially in¬ 
dependently for rather small separations, |R| > 4y2ao. 
This is in accord with recent theoretical work,^^ as 

well as experimental observation,^’^ 1 , 21,24 systems 

with small average impurity separation appear to be de¬ 
scribed in terms of the independent impurity paradigm. 
From a real-space perspective this might seem puzzling. 
A so-called ‘Kondo cloud’ of conduction electrons sur¬ 
rounds each impurity and is responsible for screening 
its magnetic moment. This screening cloud is typically 
large, characterized^^by a length scale ^ 

hvF (with vf the Fermi velocity). Our results 
show that such Kondo clouds can be substantially in¬ 
terpenetrating, even though the screening itself remains 
independent.^^ 


C. Effect of impurity hybridization strength 

As noted in Sec. HA, the impurity-host hybridization 
(and hence single-impurity Kondo temperature) can vary 
widely in real systems. Here we examine how the physics 
described in the previous subsection evolves with U/Tq. 

In Fig. 2, we consider again the local impurity spec¬ 
trum —FoImGy(ci;) at T = 0 with impurity separation 
Rcc = Ry = l,2,3,4ao, for the range of U/Tq given in 
Table I (keeping U = 3t fixed). With increasing hy¬ 
bridization To (and hence increasing we find that 

independent Kondo physics becomes relatively more im¬ 
portant, overwhelming the effects of inter-impurity in¬ 
teractions at progressively smaller impurity separations 
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R/ao = (0,1) 


R/ao = (0, 2) 


R/ao = (1,1) 



FIG. 3. QPI for two equivalent mag¬ 
netic impurities on the 3d cubic (100) 
surface, separated by vectors R/ao = 
(0,1), (0, 2) and (1,1), as shown in the 
left, center and right column panels, 
respectively. The full QPI |Ap(q, ct;)| 
is shown in the upper row panels, 
while the local and nonlocal contri¬ 
butions, obtained from Eqs. 22, 23, 
are shown in the center and lower 
row panels. Obtained for U/Fq = 28 
as per Fig. 1, with scanning-energy 
a; = 10-H = O.IT^^^P. 


(cf Eqs. 18, 19). For U/Fq = 12, the impurities be¬ 
have essentially independently already when separated 
by R = (2ao,2ao). 

The interplay between Kondo and RKKY physics is 
therefore more pronounced for larger U/Fq where the 
single-impurity Kondo scale is smaller. 


D. Two-impurity QPI 


Scattering from two magnetic impurities produces rich 
structure in the QPI. Two-impurity QPI differs from the 
idealized single-impurity limit in three key ways. First, 
the existence of a real-space vector R = r 2 — ri con¬ 
necting impurity positions gives rise to additional com¬ 
plex phase factors in the surface scattering t-matrix, see 
Eq. 12. Only in the single impurity case can these be 
eliminated by a suitable basis choice. Second, there are 
both local {a = /3) and non-local {a ^ p) contributions 
to the QPI in Eqs. 14, 15 due to the matrix structure of 
the impurity Green functions Gd{oj). Third, these Green 
functions themselves develop rich dynamics due to the in¬ 
terplay between Kondo and effective RKKY physics. As 
highlighted in Ref. 12, the resulting non-trivial scanning- 
energy and temperature dependence of the QPI is the 
hallmark of scattering from quantum impurities. 

For a generic system containing two equivalent impu¬ 
rities, the QPI can be decomposed into local and non¬ 
local contributions, Ap(q, cj) = A/)ioc(q, cj) + A/)ni(q, ^)- 


Eqs. 12-15 imply, 

Ap^(q,w) = -^[Q 7 (q,w) - Q^(-q,w)*] , (22a) 

(q, w) = (co) X (q, w) , (22b) 

with 7 =‘loc’ or ‘nl’ [here G;f'=(w) = Gy(w) = Gf{uj) 
and G’f{ijo) = G^{u) = G^^(a;)]. The host functions are 
defined as, 

Aioc(q, w) = cos(|q • R) X A(q, u) , 

A„i(q,w)= [ cos(iq-R-k||-R) 

J 

IBZ 

xG°(k||,w)G°(k|| -q,a;) , 


(23a) 

(23b) 


where A(q, ui) appearing in Eq. 23a is the single-impurity 
result, 

A(q,‘^)= / <^0(^11 >‘^)G'o(k||-q,w). (24) 

J ^^BZ 

IBZ 

Interference between scattering events from each impu¬ 
rity leads to a modulation of the QPI, controlled by the 
inter-impurity separation vector R. In particular, the lo¬ 
cal contribution Apioc(q, cj) in Eq. 23a carries an overall 
factor cos(^q-R) which lowers the symmetry of the QPI, 
and introduces nodal lines of suppressed scattering when 
q • R = (2n + l)7r. This is demonstrated in the center 
row panels of Fig. 3, where color plots of |Apioc(q)| are 
presented for R/ao = (0,1), (0,2) and (1,1) in the left, 
center and right columns, respectively. 
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By contrast, the relative scattering phase does not ap¬ 
pear simply as an overall prefactor to the non-local con¬ 
tribution Ap^\{q^,uj) in Eq. 23b. The nodal structure is 
therefore more complicated, as shown in the lower row 
panels of Fig. 3. 

Importantly, the local and non-local contributions to 
the total QPI Ap(q), plotted in the upper row columns, 
are weighted according to the local and non-local Green 
function elements — see Eq. 22. As discussed in 
Sec. IV B, the impurity separation has a major effect 
on the impurity dynamics. For R/ao = (0,1), the non¬ 
local contribution to the QPI dominates at low ener¬ 
gies uj = reflecting the strongly antiferromagnetic 

RKKY interaction which binds the impurities together 
in a singlet state. This is shown in the left column pan¬ 
els of Fig. 3 (note the rescaled color range); the full QPI 
Ap(q) therefore resembles closely its non-local contribu¬ 
tion, Apnl(q). 

For R/ao = (0,2) by contrast, the local contribu¬ 
tion dominates (again, note the rescaled color range for 
the center column panels), and the structure of the full 
QPI Ap(q) now resembles closely Apioc(q). As shown 
in the right column panels, both local and non-local ef¬ 
fects are important when the impurities are separated 
by R/ao = (1,1): One to the fine balance of Kondo and 
RKKY physics in this system. 

When the two impurities are well-separated, each is 
Kondo-screened essentially independently of the other 
(see Sec. IV B). This is reflected in the vanishing of non¬ 
local impurity Green functions 0, and hence 

the vanishing of Apni(q)- As shown in Sec. IV B this 
happens already for |R| > 4\/2ao. 


V. MANY IMPURITIES 

Solutions of full Y-impurity problems are inherently 
formidable because they involve N coupled screening 
channels. In particular, NRG (usually considered the 
numerical method of choice for quantum impurity prob¬ 
lems) cannot traditionally be used for Y > 2 impu- 
rity/channel systems. Recent developments with NRG 
— which either utilize model symmetries fully, 
or exploit the fine-grained RG description obtained by 
interleaving different Wilson chains^^’^^ — significantly 
reduce the computational cost of solving multichannel 
problems. Although models with Y = 3,4 and possibly 
5 channels can now be treated with NRG, true many- 
impurity systems remain out of reach. Furthermore, such 
real-space problems cannot generally be cast in the re¬ 
quired form involving Y impurities connected to the end 
of Y deeoupled Wilson chains. The matrix Dyson equa¬ 
tion, Eq. 4, can only be diagonalized by a single canon¬ 
ical transformation of operators in special cases where 
elements of the hybridization matrix possess the sym¬ 
metry Tii(uj) = r^°‘^(co’) and Ti^j{uj) = r^^(co’). Only 
in this case can one obtain a representation in which 
G^(uj) = (cj -h iO^ — — rjy(ci;) — Fjy)“^ are indepen- 



FIG. 4. Average fraction of impurity clusters of size Yc = 
1, 2, 3 and > 4, vs impurity coverage on a cubic lattice surface. 


dent, as required for treatment with NRG. 

However, Secs. IV B and IV C revealed that two impu¬ 
rities behave essentially independently when |R| > Rdii; 
and that impurities with realistic parameters on the pro¬ 
totypical 3d cubic lattice can be treated independently 
for a rather small separation R^ii ~ 4\/2ao. Depending 
on the specific real-space distribution of impurities, the 
system could therefore be decomposed into essentially de¬ 
coupled and hence independent smaller clusters. Indeed, 
if all impurities in a given finite sample region are sep¬ 
arated by more than Rdii, each impurity can be treated 
independently. This situation corresponds to the ‘dilute 
limit’: here an accurate description of the local physics 
of each impurity is obtained from a model involving a 
single impurity in the same host. 

For given impurity surface coverage, (j) = N/Ng (with 
Ng the number of surface lattice sites), the average inter¬ 
impurity separation (R) ^ ao(0)“^/^. But even for 
low coverage where (R) ^ Rdii, the sample may oc¬ 
casionally contain groups of a few impurities located 
close enough that their mutual interactions cannot be 
neglected. These rare ‘clusters’ may nonetheless still be 
well-separated from other impurities or clusters, whence 
independent clusters can be treated separately. At a 
given coverage, and with a specific impurity distribution, 
one can then describe accurately the many-impurity sys¬ 
tem in terms of independent impurities and clusters. 

The key question then is: what is the typical occur¬ 
rence of impurity clusters of size Yc at a given impurity 
coverage? For given coverage, and surface sample size 
relevant to FT-STS experiments, the average number of 
clusters of size Yc = 2, 3,4,... can then be determined. 

To answer this. Fig. 4 shows the configurationally- 
averaged fraction of impurity clusters as a function of 
surface coverage (with the impurities distributed at ran¬ 
dom). Specifically, we take a cluster of size Yc to com¬ 
prise Yc impurities, each within a distance Rdii of an¬ 
other element of that cluster (but separated from other 
impurities by more than Rdii). In the following we use 
Rdii = 4\/2ao, although we note from Sec. IV C that this 
is a relatively conservative definition of the cluster: larger 
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leads to smaller R^ii, and hence the earlier onset 
of the dilute limit. 

At a representative impurity coverage of 0.04%, we find 
that ^ 96.9% of impurities are independent (i.e. a ‘clus¬ 
ter’ size Nc = 1), 3% must be treated as an A'c = 2 
impurity pair (c.f. Sec. IV), and 0.1% are in clusters of 
Nc > 3. In a 500 x 500 sample region, 0.04% impurity 
density corresponds to 100 surface impurities: the surface 
is already rather crowded, with (R) ^ 28ao. But even 
then, there are on average only 3 clusters with Vc = 2 in 
such a sample. 

We note that impurity coverage can be controlled 
through various experimental protocols. For example, 
the extreme dilute limit of 10“^% impurity coverage can 
be accessed (see e.g. Ref. 54), and even specific distri¬ 
butions can be realized by manipulation of individual 
impurities.Doped thin films (such as Cr on Cu in 
Ref. 55) realize effective 2d systems with impurity densi¬ 
ties on the order of ^ 0.01%. Alternatively, evaporation 
techniques from wires onto surfaces (see e.g. Refs. 3 and 
11) controllably achieve up to ^ 0.1% impurity density; 
while dilute alloys^^ are often characterized by an equiv¬ 
alent surface impurity coverage of 0.03 — 0.3%. Higher 
impurity coverages are more common for certain transi¬ 
tion metal surfaces, which can be difficult to clean. 

The results of Fig. 4 show that at these experimentally- 
relevant surface coverages, there is a clear scale separa¬ 
tion for cluster occurrence. This suggests a hierarchy of 
approximation in which the contribution from very rare 
clusters with size Nc > N^^^ is neglected. Interestingly, 
at the above typical coverages, many-impurity systems 
are very well described in terms of independent impu¬ 
rities and independent clusters of size Nc = 2. Such 
systems are therefore amenable to treatment with NRG, 
using only single-impurity results (as in Ref. 12) and two- 
impurity results (as in Sec. IV). 

While the above analysis assumes a random impurity 
distribution, we add that inter-absorbate interactions in 
some systems may favor additional clustering. Although 
the impurity dimer or trimer occurrence for a given av¬ 
erage % coverage might then be enhanced over the pre¬ 
diction of Fig. 4, the independent cluster approximation 
nevertheless remains valid. 


A. Independent impurities: dilute limit 





FIG. 5. Independent magnetic impurities deposited randomly 
on a 3d cubic lattice surface. Upper panels: real-space LDOS 
maps. Lower panels: corresponding QPI. Left panels for 
N = 100 impurities within a 500 x 500 sample region, while 
right panels show a 100 x 100 sub-block containing V = 5 
impurities. Calculated for U/Fq = 28 and scanning-energy 
a; = 10-% = O.IT^'^^U 


physics at this level. Importantly, off-diagonal elements 
of the t-matrix thus vanish, Tq,^(cj) = = 

cj), and the QPI comprises only ‘local’ con¬ 
tributions (c.f. Eqs. 22, 23 for the two-impurity case). 

The impurity distribution does however strongly af¬ 
fect the real-space LDOS map through Eq. 6, even for 
independent impurities. The resultant QPI thus also de¬ 
pends on the impurity distribution. It is given generally 
by Eq. 14, but with (c.f. Eq. 15) 

Q(q, w) = T(q, w) x A(q, w) (25) 

where A(q, uj) is the single-impurity function Eq. 24; and 
- reflecting the independent and identical nature of the 
impurities - the full t-matrix T(q, cj) is the product of 
a structure factor, S{q) = local t- 

matrix Tioc(<^) = 

T{q,uj)=S{q)xT,M (26) 


Motivated by the above, we first consider a many- 
impurity system within an independent impurity picture. 
Specifically, we study the semi-infinite 3d cubic lattice, 
with N = 100 magnetic impurities placed randomly in a 
500 X 500 sample plaquette on the surface (corresponding 
to 0.04% impurity density). 

The independent impurity approximation implies 
a strictly diagonal hybridization matrix r(ci;) = 
U^G^(ri, r^, cj)!. Each impurity then behaves identically 
to a single quantum impurity in the same host: the 
impurity distribution does not affect the local impurity 


The experimentally-measurable power spectrum of the 
QPI follows as 

|A/)(q,w)|2 = (y;cos[(r„ - • q] ^ x|A/)iimp(q,a;)|^ , 

a,l3 

(27) 

where A/9iiinp(q, w) = [G'Jf'=(w) x A(q,a;)] is the 

pristine single-impurity QPI.^^ As such, the QPI for a 
many-impurity system in the dilute limit is related simply 
to the single-impurity QPI, but with an overlaid moire 
pattern due to the superposed cosine factors in Eq. 27. 


|Ap(q)|/ViV |Ap(r)| 
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FIG. 6. Left panel: block-averaged and symmetrized QPI, 
obtained by combining M = 25 sub-blocks of size 100 x 100 
from the original 500 x 500 sample shown in the left panels 
of Fig. 5. The right panel shows the pristine single-impurity 
result for comparison. 


I^Plimp(q) I 





This is demonstrated in Fig. 5, for an A/" = 100 impu¬ 
rity system. The upper left panel shows the full spatial 
map of the LDOS, which allows the impurity positions to 
be clearly identified. Despite the comparatively high sur¬ 
face density of 0.04% and the apparent crowding, there 
are only two Nc = 2 impurity clusters in the 500 x 500 
surface sample (although recall that in this subsection 
all impurities are treated independently). The upper 
right panel shows an expanded view of the LDOS for 
a 100 X 100 sub-block sample containing N = b impuri¬ 
ties. The corresponding QPI for both samples is shown 
in the lower panels of Fig. 5. 

Although the basic structure^^ of the pristine single¬ 
impurity QPI survives, the many impurity QPI appears 
noisy due to the impurity disorder. The q-space reso¬ 
lution of the QPI in the lower right panel is naturally 
more coarse than that in the lower left, because the real- 
space sample is smaller. The moire structure of the QPI 
due to the N = b impurities is also much more apparent, 
with well-defined interference fringes characteristic of the 
underlying impurity separation vectors. 

For random independent impurities, the QPI calcu¬ 
lated via Eq. 27 from an L x L surface sample in the limit 
L ^ oc recovers the single-impurity result.Equiva¬ 
lently, the single-impurity QPI can in principle be recov¬ 
ered by averaging the QPI over different disorder realiza¬ 
tions, since 


|Apiimp(q,w)|^ 


(|Ap(q,a;)|^) 
N + Ar2(5q,o 


(28) 


where we have used (|5'(q)p) = N Such config¬ 
urational averaging may not however be experimentally 
feasible, due to the difficulty in controlling surface topog¬ 
raphy across disparate real-space regions. 

Instead, one can consider a single L x L sample region 
decomposed into M smaller sub-blocks of size x Lb. 
as in Eig. 5. The QPI from each block can be calcu¬ 
lated and the results averaged. If the underlying lattice 
symmetry is known, the QPI pattern can also be sym¬ 
metrized by averaging over equivalent scattering vectors 


|A/}(r) - Apdii(r)| |Ap(q)| - |Apdii(q)l| 



FIG. 7. A = 100 impurities within a 500 x 500 surface region 
of a 3d cubic lattice. Impurity distribution and parameters 
as in Fig. 5. Left panel: spatially-resolved LDOS difference 
|Ap(r) — Apdii(r)|, where Ap(r) is calculated within the in¬ 
dependent cluster approximation (treating inter-impurity in¬ 
teractions in the two Ac = 2 clusters exactly), while Apdii(r) 
is calculated assuming the dilute limit within an independent 
impurity picture. Right panel: corresponding QPI difference 
|Ap(q) - Apdii(q)|. 


q. This method reduces the effect of impurity disorder in 
the QPI at the expense of q-space resolution. This proto¬ 
col was used to obtain the QPI in the left panel of Fig. 6, 
averaging over M = 25 non-overlapping sub-blocks of 
100 X 100 lattice sites. The result is clearly seen to be 
quantitatively comparable to the ‘ideal’ single-impurity 
QPI, shown in the right panel. 


B. Independent clusters: beyond the dilute limit 

We now go beyond the dilute limit by taking into ac¬ 
count explicitly the contribution from impurity clusters 
whose mutual interactions are important. These clusters 
are nevertheless well-separated from other impurities or 
clusters, and therefore remain independent. 

Within this approximation, the hybridization matrix is 
5/ocA;-diagonal, with taken to be non-zero only 

for impurities a and (3 that are both members of the 
same cluster. The t-matrix Ta^{uj) then has the same 
block diagonal structure, and the full QPI is as usual 
computed using Eqs. 14 and 15. 

As a concrete example, we use the same impurity dis¬ 
order realization as in Fig. 5, corresponding to A = 100 
magnetic impurities within a 500 x 500 region on the 3d 
cubic lattice surface (the specific distribution can be vi¬ 
sualized through the LDOS modulations in the upper left 
panel). In this representative case there are two Ac = 2 
clusters, and no larger clusters. The two-impurity clus¬ 
ters are treated independently from the remaining im¬ 
purities, which are themselves taken to be independent 
single-impurity problems. Within this framework, the 
quantum impurity problems are solved using NRG, and 
the QPI calculated. Since inter-cluster correlations are 
negligible when clusters are well-separated (Sec. IV B), 
this independent cluster approximation is essentially ex¬ 
act. 
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The local impurity physics of the Nc = 2 clusters is 
strougly affected by iuter-impurity iuteractious (Sec. IV), 
compared with isolated or iudepeudeut impurities. This 
effect would be observable through spectroscopic local 
measuremeuts with STM, albeit that the coutributiou to 
the overall QPI is expected to be weak because scatteriug 
is domiuated by the 96% of the impurities iu the sample 
which are effectively iudepeudeut. 

The left pauel of Fig. 7 shows the difference betweeu 
the real-space LDOS calculated usiug the iudepeudeut 
cluster method, aud the LDOS obtaiued through the iu¬ 
depeudeut impurity approximatiou of Sec. V A. We hud 
that this differeuce is highly localized to the cluster sites 
themselves (uote the logarithmic color scale). This cou- 
firms that bulk properties are uot strougly affected by 
the coutributiou from dilute clusters (although au accu¬ 
rate descriptiou of local quautities does of course require 
the more sophisticated treatmeut of iutra-cluster iuterac¬ 
tious). 

Iu the right pauel of Fig. 7 we plot the correspoudiug 
differeuce iu the QPI, as calculated usiug the iudepeu¬ 
deut cluster aud iudepeudeut impurity approximatious. 
The chauge iu the QPI is domiuated by iutra-cluster cou- 
tributious for au Vc = 2 pair with separatiou vector 
R/ao = (2,2), as seeu directly from the modulatiou of 
QPI iuteusity iu the Figure (the impurities iu the other 
Nc = 2 cluster are separated by R/ao = (4,4); iutra- 
cluster iuteractious are weak iu this case, aud so its cou¬ 
tributiou to the overall QPI is small). 

Our results iudicate that the iudepeudeut impurity pic¬ 
ture, strictly applicable ouly iu the dilute limit, works 
surprisiugly well wheu the cluster iucideuce is low. Theu 
the physics of iudepeudeut quautum impurities domi- 
uates the QPI — eveu for relatively large impurity deu- 
sities up to 0.1% wheu sampliug a surface regiou of size 
L X L, with L ^ 100 — 1000. To capture iuter-impurity 
iuteractious, which become importaut for local quautities 
uear clusters aud for hue details iu the QPI, the iudepeu¬ 
deut cluster method cau be employed. 

VI. CONCLUSION 

We have studied systems of multiple maguetic impuri¬ 
ties deposited ou the surface of a 3d metallic host, aud the 


associated scatteriug siguatures iu QPI. From solutiou of 
the real-space two-impurity model, we are able to defiue 
aud ideutify a ‘dilute impurity limit’ for mauy-impurity 
systems, iu which each impurity iu practice behaves iu- 
depeudeutly of the others. For realistic parameters, the 
leugth scale Rdii for the ouset of the dilute limit is fouud 
to be strikiugly low — the iuter-impurity separatiou ueed 
ouly be a few lattice sites. Overall, the physics is over- 
whelmiugly domiuated by iudepeudeut siugle-impurity 
effects for impurity coverages up to ^ 0.1% wheu the 
Koudo temperature is small, Tx ^ O.IK; while for 
Tk ^ lOOLC, a surface coverage up to ^ 1% is fouud 
to remaiu ‘dilute’. Measurable QPI therefore reflects 
the pristiue siugle-impurity result,modulated ouly by a 
trivial structure factor due to the real-space impurity dis- 
tributiou. As such, the temperature aud scauuiug-euergy 
depeudeuce of the QPI are eutirely characteristic of the 
uuderlyiug siugle-impurity Koudo effect. 

Goiug beyoud the dilute limit to higher impurity cover¬ 
age, iuter-impurity iuteractious withiu impurity clusters 
must uaturally be takeu iuto accouut. The clusters them¬ 
selves cau uevertheless be well-separated, aud therefore 
behave iudepeudeutly. Au ‘iudepeudeut cluster’ approx¬ 
imatiou, applyiug iu a ‘dilute cluster limit’ theu yields 
accurate results for the physics of such mauy-impurity 
systems. 

Au iudepeudeut cluster picture must of course ulti¬ 
mately fail wheu impurity coverage is iucreased towards 
aud above a percolatiou threshold where the eutire sys¬ 
tem becomes ‘couuected’. The physics theu iuvolves a 
subtle iuterplay betweeu Koudo, RKKY, disorder, aud 
the collective heavy fermiou physics of the diluted peri¬ 
odic Audersou lattice - aud remaius a pereuially opeu 
problem. 
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